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Abstract 



^ , We study the properties of the vortex state of a trapped Bose gas in the large-gas-parameter 
> regime. To test validity of the Gross-Pitaevskii theory in this regime for the vortex states we 

o ; 

compare the results of the Gross-Pitaevskii and the modified Gross-Pitaevskii calculations for the 
I W e„e.„ the *e.n.a, po.e.i., .he de„... p^ohle ... .he shift o. .he ,.a<i.pole 

I 

modes of the collective oscillations of the condensate. We find that in the large-gas-parameter 
Q regime two calculations give substantially different results for all the properties mentioned above. 

^ PACS numbers: 03.75.Fi, 03.65.Db, 05.30.Jp 
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I. INTRODUCTION 



The mean-field Gross-Pitaevskii (GP) equation has been quite successful in describing 
both static and dynamic properties of Bose-Einstein condensates (BEG) of alkali atoms 
confined in magnetic or optical traps 

mm 

. It is well known that the GP theory is 
valid for condensates satisfying dilute gas condition na^ << 1 ( where n is atomic density, a 
is the s-wave scattering length of interatomic potential and the parameter x = na^ is called 
the gas parameter). Physically it implies that an interacting atomic gas can be considered 
to be dilute when the average distance between the atoms (oc n~^/'^) is much larger than the 
range of the atomic potential (oc a). Moreover,under this condition the detail shape of the 
inter-atomic potential is not required rather the boson-boson interaction can be simulated by 
a pseudo potential which is completely characterised by scattering length a. It is important 
to note here that in some recent theoretical studies it has been shown that as long as 
X < 10~^ the GP theory produces accurate results for both static B y, y, 0|and dynamic 



properties 
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12l |. Typically the value of the gas parameters in most of the experiments 
were in the range 10~^ — 10~^ and therefore the GP theory works very well in predicting the 
properties of these trapped condensates. However, in a recent experiment condensates with 
peak gas parameter of the order of 10~^ have been achieved by enhancing the scattering 
length a with the help of Fesbach resonance [isl. For condensates with such large values of 



the gas parameter the question of validity of GP theory has been raised and tested [10|, [12 1 
by employing modified GP (MGP) theory within the local-density approximation. These 
studies clearly show that for large values of gas parameter 10^^) both ground-state 
properties ( energy, chemical potential) and djTiamic properties ( frequencies of collective 
oscillations) are significantly modified by MGP theory. 

The vortex states in trapped condensates play an important role in establishing the 
superfiuid properties of the condensates . It is then natural to ask how the properties 
of vortex states are modified in the large-gas-parameter regime. This has motivated us 
to study the properties of the trapped condensates with vortex states in the large-gas- 
parameter regime. To describe the condensate in the large-gas-parameter regime we employ 
the MGP theory which is considered to be applicable in this regime. By comparing the 
MGP results with the corresponding GP numbers we make a systematic assessment of the 
corrections introduced by the MGP equation over the GP results. In this paper we consider 
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condensates with the values of gas parameter spanning a range similar to that achieved in 
the experiment of Cornish et al. |l3|. The results presented in this paper can be divided 
in to two parts. In the first part we focus our attention on the modifications of energy, 
chemical potential and the density profile of the condensate with a quantized vortex. 

The quantized vortex states not only affect the static properties of the condensate but 
the dynamic properties also get modified by them. For example, the presence of a quan- 
tized vortex state leads to the splitting of the frequencies of the two modes of quadrupole 
oscillations with opposite values of the third component of the angular momentum, which 



le angular r 
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are degenerate in the absence of a quantized vortex state 
in a condensate breaks the time reversal symmetry which in turn results in the removal 
of the degeneracy of the two modes of oscillations carrying opposite values of angular mo- 
mentum. Therefore, the splitting of the two quadrupole modes of collective oscillations can 
be employed to detect the presence of a quantized vortex states in the BEC 1^1 as the 
measurement of the frequencies of the collective oscillations can be carried out with high 
precision . In the second part of this paper we focus our attention on the effect of the large 
gas parameter on the splitting of two degenerate quadrupole modes due to the presence of 
a vortex state. We calculate and the compare the results for the frequency shift of the two 
quadrupole modes by employing the GP and the MGP equations in the large-gas-parameter 
regime. 

For the purpose of calculation we use the variational approach with quite accurate vari- 
ational ansatz for the wave function of a trapped condensate with a quantized vortex line 
along the z-axis and the frequency shifts of the quadrupole modes are then obtained by 
employing well estabhshed sum-rule approach jl^ ll^ • 

The paper is organized in the following manner. In section II we describe the theoretical 
methods employed in this paper and briefly describe the MGP theory followed by the vari- 
ational method employed to obtain the wave function of the vortex state and other physical 
observables mentioned above and the sum-rule approach for the calculation of frequencies of 
collective oscillation. The section III is devoted to the discussion of the results. The paper 
is concluded in section IV. 
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II. THEORY 



The ground-state energy functional associated with a condensate of N bosons each with 
mass m confined in a trap potential Vt{r) can be written as 3] 



dr 



;^|V^P + V;(r)|^|2 + eH|$| 
2m 



(1) 



where \l/(r) is the the condensate wave function (order parameter) and n{r) represents the 
corresponding density and it is given by n(r) = |^(r)p . The condensate wave function 
\E'(r) can be determined by minimizing the above energy functional. In the above equation 
the first, the second and the third terms represent the kinetic energy of bosons, the energy 
due to the trapping potential, and the inter-atomic interaction energy within the local den- 
sity approximation (LDA), respectively. To go beyond the GP theory we make use of the 
perturbative expansion for e{n) in terms of the gas parameter na^, 
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(2) 



The first term in the above expansion, which corresponds to the energy of the homogenous 
Bose gas within the mean-field theory as considered in the GP theory, was calculated by 
Bogoliubov i21j. The second term was obtained by Lee, Huang and Yang (LHY) 22], while 
the third term was first calculated by Wu [2^ using the hard-sphere model for the interatomic 
potential. Although, it has been emphasized in the literature that the above expansion is 
valid only for na^ << 1, it is only recently that the range of validity of above expansion has 
been systematically investigated by Giorgini et al. 24|. They have used diffusion Monte 
Carlo (DMC) method to calculate the ground-state of uniform gas of bosons interacting 
through different model potentials with a > 0. It has been found that the expansion (j21) 
is valid as long as na^ < 10^'^. However, for gas parameter beyond 10"'^ the inclusion of 
the logarithmic term leads to a severe mismatch with the DMC simulation results. On 
the other hand, expansion (0) up to the LHY term gives an accurate representation of the 
DMC calculations even for the gas parameter of the order of 10~^. Consequently, we do not 
consider the logarithmic term in the expansion (j2)) for all the calculations in this study. 

The trapping potential Vt{r) is taken to be axially symmetric characterized by two angular 
frequencies tu^ and u'^ {u^ = cu^ = cu^ 7^ c<j°). It is given by 

,o2 



VAr) 



muj , 



(3) 



where Aq = uj^^/uj\_ is the anisotropy parameter of the trapping potential (Aq = 1 corresponds 
to a spherically symmetric trap). 

The ground-state energy functional (Q) can be easily generalized to include the vortex 
states. In this paper we consider condensates with a vortex line along the z-axis and all the 
atoms are flowing around it. The wave function for such a condensate can be written as |25] 



^(r) = V(r)e*^'^ 



(4) 



where (f) is the angle around the z-axis and k is an integer denoting the quantum circulation 
and the total angular momentum along the z-axis is given by Nhn. Substituting the above 
complex wave function of Eq.(j3]) in Eq. ^ we get generalized energy functional for the 
condensate with a vortex state as Ul 



E\ 



dr 



2m 2mr , 



(5) 



In the above equation r± = \/x^~-\-~y^ and the density n{r) = |'?/'(r)p. The the presence of 
centrifugal term due to the vortex state makes the above functional different from Eq. (^. 
The minimization of the above functional with respect to ilj{r) with the constraint 



leads to the MGP equation for the condensate with vortex state 

32a3/2, 



(6) 



V H ' ^^^-^ ' 
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(7) 



where /i is the chemical potential arising from the constraint condition given by Eq. The 
GP equation for the vortex state can be obtained from Eq. (0) by neglecting the interaction 
energy term proportional to in this equation. It is not possible to find an exact solution 
to Eq. (jZj). Consequently various numerical techniques have been developed to solve the 
above nonlinear Schrodinger equation for studying the properties of trapped BEG j^^El- 
In this paper we use the variational method to obtain the wave function (or density) and 
other relevant physical observables of the condensate. The main advantage of this method 
is that with a suitable choice for the variational form for the wave function one can get 
,.ite acc.a. ..1. wi.h le. co.pu.a^onal effo.. In eaHie. wo.. MU^^ have 
demonstrated the applicability of the variational method by calculating the ground state 
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properties and the collective oscillation frequencies of condensates both with the GP and 
the MGP equations for wide range of values of particle number and scattering length. 

Following Ref. the variational form for the wave function of the condensate with 
quantized vortex state is chosen as 

iI}{yi) = Arl^e V^x/ , (8) 

where g, A, uo^ and p are the variational parameters which are obtained by minimizing the 
energy functional given by El with respect to these parameters. In the above equation we use 
the scaled co-ordinate ri = r/aho with aho = (h/muj^)^. We note that n(ri) is normalized 
to unity and from Eq. (jH)) it can be seen that 

ra(r) = — n(ri). (9) 

The above variational form of the wave function has already been shown to describe the 
state of the dilute Bose gas with a q uantized vortex confined in a trap quite accurately for 
a wide range of particle numbers 20| . In this paper we extend the applicability of this wave 
function in the large-gas-parameter regime also. 

The proportionality factor A in Eq. (jH)) is determined by the normalization condition 
(Eq. ©) 



yAr(3/2 + g) /^y='/'+^) 



^ 7r3/2r(l + g)r((3/2 + q)/2p) U° / ^^^^ 
Using this variational form for the wave function physical observables can be expressed 
analytically in terms of the four variational parameters. For example, the energy functional 
El = E/huj"^ can be written as 

^ = T + U + Erot + El, + El,. (11) 

Here T, U and Erot denote the average kinetic, the trapping potential and the rotational 
energies per particle, respectively. The term E^^^j gives the interaction energy per particle in 
the mean-field approximation as considered in the GP theory while Ef^^ gives the correction 
due to the LHY term in the expansion Q resulting in the MGP equation. The analytical 
expressions for these energy components are 

+ 2q) [(1 + 2p)(l + A/2) + q{2p + 2g + A)] r((l + 2g)/2p) 

u;0 4(3 + 2g)r((3 + 2g)/2p) ' ^ ^ 
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ujl{l + q + \l/2\)T{{b + 2q)/2p) 

o;^(3 + 2g)r((3 + 2g)/2p) ' ^ ^ 

_ ^V(l + 2g)r((l + 2g)/2p) 

4a;0gr((3 + 2g)/2p) ^'^^ 



pi _g..- pv/Ar^(g + 3/2)r(2g + l)r((4g + 3)/2p) 

2(4'?+3)/2p0Fr2((2g + 3)/2p)r2(g+l)r(2g + 3/2)' 



(15) 



and 



9/4 /o\ (5g+3/2p) 

p3/2;^3/4 



p2 _ 128 2v/2 / /,/^^y/Y2\ 

r5/2(g + 3/2)r(5g/2 + l)r((5g + 3)/2p) 



r5/2((2g + 3)/2p)r5/2(g+l)r((5g + 3)/2) ^^^^ 
where a = a/aho and r(m) is the gamma function. It is easy to check that above energy 
components correctly reduce to the expressions derived in |9| for the case of condensate 
without quantized vortex by setting k and q equal to zero. For particular values of and d 
the parameters uj±, q, A and p are determined by minimizing the energies given by Eq. (fTT|l 
and Eqs. (fT ^ - (fT^ . 

Having described the variational method for obtaining the wave function and the other 
physical properties we now describe the sum-rule approach of many-body response theory 
for the calculation of frequencies of collective oscillations. In the coUisionless regime the 
collective excitation frequencies of a confined bosonic gas is well described by the sum 
rule approach. The most important advantage of this approach is that the calculation 
of frequencies requires the knowledge of the ground-state wave function (or the ground- 
state density) of the many-body system only. This method has been extensively applied to 
calculate the frequencies of the collective oscillations of trapped atomic gases jOj. In this 
paper we follow the sum-rule approach developed in jlj] for the calculation of the splitting 
of the quadrupole modes of collective oscillations. 

The collective excitations of any many-body system are generally probed by applying 
an external excitation field. For excitation operators and F_ exciting two quadrupole 
modes with o ppo site angular momentum, respectively, the corresponding strength functions 



are given by 



IS 



S^{E) = Y,mF±\n)\'6{E - E^) (17) 
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where |0) denotes the ground state of the system and En — Eq is the excitation energy of 
the eigenstate \n) of the Hamiltonian H relative to the ground state energy Eq. The p-th 
order moments of the strength function can be defined as 

= j EP{S+{E) - S4E))dE. (18) 

These moments provide various energy weighted sum rules which are employed to obtain 
the frequencies of collective oscillations. One of the important properties of these moments 
is that, for a given p, some of the moments can be expressed in terms of the commutators of 
the excitation operators with the many-body Hamiltonian H. Some of the energy weighted 
sum rules which are relevant for this paper are 

mo = (0|[Ft,F]|0) 

mt = (0| [F\[H,F]] |0) 

= {0\[[F\H],[H,F]] |0) 

mt = {0\[[F\H],[H,[H,F]]]\0) (19) 

where |0) represents the ground state of the hamiltonian H and [, ] represents the commutator 
between the corresponding operators. 

These expressions can be used to calculate the frequency shifts of the quadrupole oscilla- 
tions of a trapped condensate. For the quadrupolar case we will consider the modes excited 
by the operators 

F± = (x±zi/)2 (20) 

carrying angular momentum = ±2, where represents the z-component of the angular 
momentum of the elementary excitation. In order to use the above expressions for the 
calculation of the shift in frequencies of quadrupole modes we now assume that the moments 
mp are exhausted by a single excitation with frequency uj± and strength a±. Under such 
single-mode approximation the strength distribution can be written as 

S±{E) = a±5{E-nuj±). (21) 

Moreover, the vanishing of mg ( as given by the first moment of Eq. (jl9|) ) due to the 
commutator of F+ and F_ being zero leads to the result (t+ = (t_ = a. We note here that 
the single mode approximation is well suited for the condensates with large number of atoms 
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N and positive scattering length a 13]. Using Eqs. (PT|) and (fTHj) we obtain the following 
expression for the difference between the two frequencies 

h6 = h(uj+-uj) = (22) 

Furthermore, the individual frequencies of the quadrupole oscillations can be determined by 
using Eq. ()22|) along with the expression for the square of the mean frequency given by 

In the next section we calculate the moments , and in terms of some averages 
over the vortex state wave function and use Eqs. ()22j) and ()23|1 to determine the frequency 
shifts of the two quadrupole modes due to the presence of a quantized vortex in a trapped 
condensate falling in the large gas parameter regime. 

III. RESULTS AND DISCUSSION 
A. Ground state properties 

We begin this section with a comparison of our results obtained by the variational method 
with the corresponding numbers of Ref. jbll calculated by solving the MGP equation by 
steepest descent methods. We make this comparison to check the accuracy of the variational 



approach as it has been shown in Ref. |3l| that their MGP results are in very good agreement 
with the results of ab-initio variational Monte Carlo calculation both in the presence and in 
the absence of vortex state. At this point we note that in addition to the comparison with 
the ab-initio results the virial relation among the different energy components also provides 
a way of checking the correctness and the accuracy of the variational solutions. For the 
MGP theory in the presence of a single quantized vortex the virial relation is given by 

2T-2U + 2Erot + ^Elt + \eI, = (24) 

This expression has been derived by using the variational nature of the ground state energy 
and the scaling transformation n(r,,) — > if'n{Y). In all our calculations the virial relation 
given above is satisfied up to sixth decimal place or better. This indicates that the the 
variational method employed in this paper yields quite accurate results for condensates with 
vortex state in the large-gas-parameter regime. 



For the purpose of comparison we calculate the total energy, individual energy compo- 
nents given by Eqs. (fT^ - IT^ and the chemical potential of a ^^Rb condensate carrying a 
single vortex state (k = 1). In accordance with the Ref. ^ we consider a disk-shaped 
trap with Aq = and cj^ = 27r x 220 Hz confining = 500 condensate atoms and the 
scattering length is chosen to be a = 0.15155aho- The results of this calculation are shown 
in Table I along with the corresponding numbers of Ref. [sij] in parentheses. It can be 
clearly seen from Table I that numbers obtained by the variational method are very close to 
the corresponding results of Ref. In particular our numbers for the total energy, trap 

energy and both the components of the interaction energies are slightly higher (of the order 
of 1%) than those of Ref. On the other hand, our results for the kineteic energy and 

the rotaional energy are appreciably lower and higer respectively, than the corresponding 
energies reported in Ref. j3l|. Despite the difference between the results for the kinetic and 
the rotational energies the results obtained for the total energy and the chemical potential 
obtained by the two metods match quite well due to cancellation of errors. The cancellation 
of errors also lead to very good satisfaction of the virial relation. Moreover, we find that 
for the above choice of and a the MGP reuslts for the total energy and the chemical 
potential are approximately 19% and 16% respectively, higher than the corresponding the 
GP numbers. These results along with our 

applicability of the variational approach in calculating the properties of condensates both 
with and without quantized vortex state in the large gas parameter regime. 

Now we focus our attention on the application of the MGP theory to condensates with 
a single quantized vortex for larger values of gas parameter as achieved in the experiment 
of Cornish et al. 3|- In accordance with the experiment we consider = 10*^ ^^Rb atoms 
confined in an anistropic trap with frequencies a;^^ = 2it x 17.5 Hz and uj^ = 2Tr x 6.9 Hz 
and the scattering length is varied from a = 1400ao to a = lOOOOao, where is the Bohr 
radius of hydrogen atom. We note here that this range of values of the scattering length 
falls in the large gas parameter regime as the maximum value of a corresponds to the peak 
gas parameter Xpeak ~ 10~^. In Table II we show the results for the total energy and the 
chemical potential obtained by employing the MGP equation and compared them with the 
corresponding numbers calculated by the GP equation. It clearly shows that the difference 
between the MGP and the GP results increases with increase in the scattering length. For 
a/ao = 1400 the differences in the total energy and the chemical potential are approximately 
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3% and 4% respectively. On the other hand, for the maximum value of a/a^ = 10,000 the 
MGP reults for the total energy and the chemical potential are higher by around 26% and 
30% respectively, over the corresponding GP numbers. From these results we conclude 
that in the large gas parameter regime the MGP results for the condensates with a single 
quantized vortex are substantially altered as compared to the corresponding GP results. It 
is important to note here that although the MGP corrections to the energy are significant, 
the critical angular velocity Qc for the formation of the vortex state given by 2^ 

% = m/N)^ - {E,/N),] , (25) 

changes nominally as compared to the GP calculations. For example, corresponding to 
a/oo = 10,000 we obtain = 0.325uji^^ and Qc = 0.335co'^^ for a single vortex with k, = 1 
by employing the MGP and the GP calculations respectively. This observation regarding 
the small change in value of critical angular velocity was also made in Ref . [sil , however, for 
much lower value of the gas parameter. Our results clearly shows that the MGP calculation 
does not provide any sizable correction to the critical angular velocity Qc even for very large 
value of the gas parameter. 

Having estimated the corrections introduced by the MGP calculation over the GP results 
for the energy and the chemical potential now we wish to present the results for the density 
profile obtained by the MGP and the GP calculations and compare them to study the effect 
of large gas parameter on the density profile. In Fig. 1 and 2 we show the density profile 
of the condensates with a single vortex for the values of scattering length a/oo = 1400 
and a/aQ = 10000 respectively obtained by the MGP and the GP calculations. It can 
be clearly seen from these two figures that for a/ao = 1400 the MGP calculation leads 
to slight modification in the density profile over the GP result. On the other hand, for 
a/ao = 10000 the MGP density profile is substantially different from the corresponding GP 
density profile. Due to the repulsive nature of the LHY term in the MGP equation the 
bosons in the condensates are pushed outward resulting in the increase in the size of the 
condensate and decrease in the value of maximum density in comparison to the GP case. 
Furthermore, it is evident from Fig. 2 that the values of the peak density npeak and the the 
transverse distance r^^"^ at which densities become maximum are significantly different in 
the GP and the MGP calculations . To show the difference in the MGP and the GP density 
profiles in more quantitative manner we plot in Fig. 3 and 4 the peak density rip^ak and the 
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corresponding peak position r^^*^ respectively, as functions of the scattering length. From 
the variational form of the wave function given by Eq. it is easy to verify that the density 
peaks at 



The values of variational parameters tu^, X, p and q obtained by minimization of the energy 
are used to calculate r^*^^^ and Upeak by using Eqs. (j26|l and (j27jl . From Fig. 3 we observe 
that with both MGP and the GP calculations the peak density shows a decreasing trend with 
increasing a/ao due to greater repulsion between the bosons. As mentioned above for the 
MGP calculation the second term in inter-atomic potential corresponds to more repulsion 
which results in higher reduction in the peak density as compared to the GP case. The 
difference between the MGP and the GP peak densities increases with increasing scattering 
length. For example, as a/aQ is scanned from 1400 to 10000 the difference between the GP 
and the MGP results grows from around 7% to 37%. In contrast to this r^^"^ increases as 
the scattering length is increased in conformity with the fact that the size of the condensate 
grows as the repulsive interaction goes up. However, like peak density the difference between 
the results for rf^^^ by the MGP and the GP calculations also increases with increasing value 
of a/aQ. For a/a^ = 1400 and a/a^ = 10000 the differences in the MGP and the GP r^^"'^ are 
found to be approximately 3% and 22% respectively. Furthermore we also notice from Fig. 
1 and 2 that the depth of the vortex hole around the z-axis obtained by the GP calculation 
is significantly greater than the MGP case in the large-gas-parameter regime. However, the 
results for the size of the hole around z-axis characterised by the ratio r\^^^ / a/ (0|x^ + ?/^|0) 
3^ . obtained by the two calculations show difference of lesser degree as the scattering length 
is increased. For example corresponding to a/ao = 1400 and a/ao = 10000 the differences 
in the size of the holes are of the order of 1% and 4.5% respectively. 

Thus we conclude that in presence of vortex the MGP results for the properties like 
the total energy, chemical potential and the total density profile except the critical angular 
velocity VLf. and the size of the vortex hole around the z-axis differ significantly in comparison 
to the GP predictions in the large-gas-parameter regime. In the next section we discuss the 




(26) 



and the expression for the peak density is given by 




(27) 
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results for the splitting in the quadrupole modes of collective oscillations due to the presence 
of vortex in the large gas parameter regime. 



B. Splitting of quadrupole modes 

To calculate the frequency shifts of the two quadrupole modes characterized by the com- 
ponents of angular momentum m^. = 2 and = —2 respectively, first we need to evaluate 
the moments mf, and . For excitation operators given by Eq. ()2()|1 the commutators 



in Eq. (jl9p can be evaluated to give following expressions for the moments IJ, |29 1 



— (0|x^ + y^|0) 
m 



m 



(0|x2 + j/2|0) + 



(28) 



where Pi is the i-th component of the linear momentum vector p and |0) represents the 
wave function of the vortex state obtained by solving the time independent MGP or GP 
equations (Eq. In deriving above equations the relation F| = F_ has been used. It 

is worth noticing that the moments up to third-order do not depend on the boson-boson 
coupling parameter (two-body interaction) explicitly. This is because the energies within the 
LDA do not contribute to moments up to third order for the excitation operators satisfying 
Now use of moments and m2 in Eqs. ()28|) and ()22|) vields the result 

2 Nhn 



S 



(29) 



where (r^) is the transverse size of the condensate determined by the average 

{rl) = J + y^) n{r)dr (30) 

It is important to note that albeit the frequency shift does not explicitly depend on the 
boson-boson coupling parameter but implicit dependence on the two-body interaction enters 
through the wave function or the density of the condensate which crucially depend on the 
nature of the boson-boson interaction. We have already observed that the transverse sizes 
of the condensate obtained by employing the MGP and the GP equations differ significantly 
in the large-gas-parameter regime. Consequently, the MGP and the GP frequency shifts of 
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the two quadrupole modes will also vary significantly in this regime. We show the result 
for the MGP and the GP frequency shifts as a function of the scattering length in Fig. 5. 
The frequency shift between the two quadrupole modes decreases with the increase in the 
scattering length due to the fact that higher values of scattering length correspond to greater 
repulsive interaction between the bosons resulting in condensates with larger values of {r\). 
For up to around a/ao = 1000 the GP and the MGP frequency shifts are nearly identical, 
however, as the scattering length increases further the difference between the two results 
also increases. For example, at the maximum value of the scattering length (a/ao = 10000) 
we find that the GP frequency shift is around 26% higher in comparison to the MGP result. 

Finally, we note that the quadrupole frequencies u!± can be obtained from Eq. and 
the expressions for the mf and given in Eq. (j28|) . For the values of the gas parameter 
considered in this paper the transverse kinetic energy (0|p^ + Py|0) is much much smaller 
than the transverse radius (0|x^ + y'^\0) of the condensate with a vortex state and thus it 
can be neglected in . As a result of this, under single mode approximation (Eq. (j2T| ) 
the expressions for the quadrupole frequencies are given by 

cj± = ^2101 ± - (31) 

From the above results we conclude that in the large-gas-parameter regime the MGP and 
the GP results for the frequency shift of the two quadrupole modes due to the presence of 
a quantized vortex and also the frequencies of the two modes are substantially different for 
the values of scattering lengths which are achieved by tuning the Fesbach resonance. 

IV. CONCLUSION 

In this paper we have studied the properties of BEG with single quantized vortex in the 
large gas parameter regime. For this purpose we have employed MGP equation which has 
been obtained by including the second term (LHY term) in the perturbative expansion of 
the inter-atomic interaction energy per particle obtained from the ground-state energy of 
uniform Bose gas. We have used the variational approach to solve the MGP equation by 
employing a suitable ansatz for the wave function of a condensate with a quantized vortex. 
The wave function obtained by variational approach is then employed to calculate the static 
properties like the total energy, chemical potential and the density profile for wide range 
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of scattering length lying well within the large gas parameter regime. The correctness and 
accuracies of our solutions are checked by verifying the generalized virial relation and also 
comparing them with the solutions obtained by the steepest descent method. By using the 
wave function of the vortex state and the sum-rule approach we have also calculated the 
frequency shift of the two quadrupole modes of the collective oscillations due to the presence 
of a quantized vortex in the large gas parameter regime. To test the accuracy of GP equation 
we have made a detailed comparison of the results of the GP and the MGP equations for 
all the observables mentioned above. We have found that in the large gas parameter regime 
MGP equation provide substantial amount of corrections in all the observables considered 
in this paper. These changes are quite sufficient to be observed in measurements involving 
condensates with large gas parameter as achieved in Ref . . The comparison of the results 
obtained in this paper with the experimental results will also test the validity of the GP 
theory for the description of such condensates. 
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TABLE I: Results for the chemical potential and the energies in units of hLiij_ obtained from the GP 
and the MGP calculations for the condensate of of = 500 '^^Rb atoms trapped in an anisotropic 
trap with Aq = \/(8) and tj^ = 27rx220ffz. The scattering length is o = a/a± = 0.15155. Numbers 
in parenthesis are results of Ref. jsil 







Ei/N 


T/N 


U/N 


ElJN 


Elt/N 


Erot/N 


GP 


13.334 


9.889 


0.3926 


5.8056 


3.445 




0.2463 




(13.187) 


(9.7836) 


(0.4251) 


(5.7427) 


(3.4039) 




(0.2119) 


MGP 


15.811 


11.443 


0.3405 


7.1248 


2.4999 


1.2454 


0.2324 




(15.623) 


(11.305) 


(0.37692) 


(7.0377) 


(2.4824) 


(1.2233) 


(0.1849) 



TABLE II: Results for the chemical potential and the total energy per unit number of bosons in 
units of hu)^ obtained from the GP and the MGP calculations for the condensate of of A^ = 10'^ 
^^Rb atoms trapped in an anisotropic trap with Aq = 0.39 and to^ = 27r x 17.5Hz. 



a/ttQ 


GP 


MGP 


/^i 


Ei/N 


Ail 


Ei/N 


1400 


10.223 


7.545 


10.621 


7.790 


3000 


13.649 


9.955 


14.898 


10.739 


8000 


19.955 


14.422 


24.798 


17.528 


10000 


21.774 


15.715 


28.223 


28.223 
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FIG. 1: Comparison of the GP (dashed line) and the MGP (sohd Une) density profiles of a vortex 
state with k = 1 as a function of the radial distance r±_ (in units of auo) from the z-axis for N = 10^ 
^^Rb atoms trapped in an anisotropic trap with Aq = 0.39 and uj'j_ = 2ttx 17.5Hz . The scattering 
length a/ao = 1400. 

FIG. 2: Same as Figure 1 but with the scattering length a/ao = 10000. 

FIG. 3: The peak density ripcak as a function of scattering length a/ao of a condensate with a 
quantized vortex of quantum of circulation k = 1. The solid line corresponds to the MGP result 
and the dashed line is obtained with GP calculation. The trap parameters are same as Fig. 1. 

FIG. 4: The peak position r^'^^ (in unit of aho) as a function of scattering length a/ao of a 
condensate with a quantized vortex of quantum of circulation k. = 1. The solid line corresponds 
to the MGP result and the dashed line is obtained with GP calculation. The trap parameters are 
same as Fig. 1. 

FIG. 5: The frequency shift (in unit of ujj_) of the quadrupole oscillations of a BEC due to the 
presence of a single quantized vortex (k = 1) as a function of scattering length a/ao. The solid 
line corresponds to the MGP result and the dashed line is obtained with GP calculation. The trap 
parameters are same as Fig. 1. 
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